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; ABSTRACT 

Context. The progenitors of supemovae type la are usually assumed to be either a single white dwarf accreting from a non-degenerate 
companion (the single degenerate channel) or the result of two merging white dwarfs (the double degenerate channel). However, no 
\^ \ consensus currently exists as to which progenitor scenario is the correct one, or whether the observed supemovae la rate is produced 

^Nj ■ by a combination of both channels. Unlike a double degenerate progenitor a single degenerate progenitor is expected to emit supersoft 

X-rays for a prolonged period of time (~ 1 Myr) as a result of the burning of accreted matter on the surface of the white dwarf. An 
argument against the single degenerate channel as a significant producer of supemovae type la has been the lack of observed supersoft 
I X-ray sources and the lower-than-expected integrated soft X-ray flux from elliptical galaxies. 

, Aims. We wish to determine if it is possible to obscure the supersoft X-ray emission from a nuclear burning white dwarf in an accreting 

single degenerate binary system. In case of obscured systems we wish to determine their general observational characteristics. 
Methods. We examine the emergent X-ray emission from a canonical supersoft X-ray system surrounded by a spherically symmetric 
S-M configuration of material, assuming a black body spectrum with Ti,i, = 50 eV and L = lO'^'erg ■ s"'. The circumbinary material 

Q , is assumed to be of solar chemical abundances, and we leave the mechanism behind the mass loss into the circumbinary region 

5_( ■ unspecified. 

Results. We find that relatively small circumstellar mass loss rates, M = 10"' - lO^'^Moyr"', at binary separations of ~ 1 AU or less, 
\ will cause significant attenuation of the X-rays from the supersoft X-ray source. Such circumstellar mass loss rates are sufficient to 

make a canonical supersoft X-ray source in typical extemal galaxies unobservable in Chandra. 

Conclusions. If steadily accreting, nuclear burning white dwarfs are canonical supersoft X-ray sources our analysis suggests that they 
can be obscured by relatively modest circumbinary mass loss rates. This may explain the discrepancy of supersoft sources compared to 
^ I the supemova la rate inferred from observations if the single degenerate progenitor scenario contributes significantly to the supemova 

, la rate. Recycled emissions from obscured systems may be visible in other wavebands than X-rays. It may also explain the lack of 

observed supersoft sources in symbiotic binary systems. 

■ 

■ Key words. (Stars:) supemovae: individual - (Stars:) binaries: close - Accretion, accretion disks - (Stars:) white dwarfs - Stars: winds, 
• ■ outflows - X-rays: binaries 

o 

1. Introduction Despite decades of intense research on the subject, the ex- 

act nature of the progenitor systems of these important astro- 
physical explosions remains unclear Carbon-oxygen WDs are 
Type la supemovae (SNe) are believed to be carbon-oxygen characteristically formed at masses much lower (~ 6 M^) than 



V . 

s_l ■ "15* A e f\\Tr\ \ \ \ r^u A lu *u * that needed for thermonuclear runaway (~ 1 .37 Mq), and there 

H , white dwarfs (WDs) close to the Chandrasekhar mass that un- . , , , , , . , , , ^, , 

. F^. ■ dergo thermonuclear runaway in their centers. The resulting ex- "° ^nown process by which an isolated sub-Chandrasekhar 

plosion produces radioactive iron-group elements, and the sub- ^^f' ^D can grow to the critical mass at which it explodes as 

. ■ f » »ui f 56ivT- u a SN la. Hence, It IS usually agreed that SNe la can only arise 

sequent decay of these, most notably of Ni, powers charac- . , . , „,tC 

teristic light curves that obey a well-known relation between !" binary systems, where a WD accretes matter from a compan- 

luminosity at maximum light and fall-off time (Phillips [l993). However, the exact method of accretion remains dis- 

A u CM T A A ^ A A- ui 1 ' putcd. Two binary progenitor scenaHos ai'e usualW considcrcd: 

As a result, SNe la are considered standardizable cosmologi- , . , , ,ot^. . , ■ , ,,,t^ 

cal candles. To this effect they have been utilized to suggest ''^^^^ degenerate (SD), m which a WD accretes mass from a 

that the expansion of the Universe is accelerating (Riess et al. """r'^^g^"^'!^^^ companion (niain sequence or giant star Whelan 

[MS Perlmutter et al. which in turn has given rise to "^J^^" ^^^^ degenerate (DD), m which two 

the paradigm of Dark Energy. Additionally, the energy release sub-Chandrasekhar mass WDs merge with a mass^or above 

of SNe la is large enough to influence the dynamics of their '° '''P^°'^' ' (WebbinkM Iben & 

host galaxies, and the nucleosynthesis taking place during the ' 

explosions is the main source of iron group elements in galactic While the DD scenario has garnered considerable attention 

chemistries. recently, the SD scenario has been the most popular scenario for 
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a long time, and much work has been done on the physics of 
this progenitor scenario (e.g. Hachisu, Nomoto & Kato ll996l l. 
It was shown by Nomoto ( 1982) that steady nuclear burning of 
hydrogen-rich material accreted from a companion onto a mas- 
sive (~ 1 Mq) WD can only take place in a fairly naiTow inter- 
val of accretion rates close to 10"^ M^yr^K At smaller or larger 
accretion rates it is unclear if the WD will be able to grow suf- 
ficiently in mass for a SN la to occur, due to possible mass loss 
from nova eruptions or stellar winds. This puts rather tight con- 
straints on the parameters of the progenitor systems. 

In the 1990's luminous supersoft X-ray sources (SSSs) were 
recognized as an important new class of X-ray source (Triimper 
et al. 119911 Greiner et al. il991i ). based on observations of the 
Large Magellanic Cloud made with the Einstein Observatory in 
the late 70's and early 80's (Helfand & Grabelskv ll981i) . Since 
then, newer generations of X-ray instruments such as ROSAT, 
BeppoSAX, XMM-Newton and Chandra X-ray Observatory have 
found similar sources in other galaxies, including the Milky 
Way. As the name suggests, SSSs are characterized by hav- 
ing a much softer spectrum than those of more commonly 
known X-ray binaries involving a neutron star or a black hole. 
Subsequently, van den Heuvel et al. (I1992I I showed that a mas- 
sive WD accreting from a companion star at the steady-burning 
rate will emit X-rays with a spectrum consistent with that ob- 
served for a certain subset of SSSs as a result of thermonuclear 
burning of the accreted material. This made SSSs interesting as 
possible SD progenitor systems of SNe la. 

If WDs undergoing steady nuclear burning on their surfaces 
look like SSSs, and if the SD progenitor scenario is the dom- 
inant contributor to the SN la rate, then we should expect - at 
least naively - to see a corresponding population of SSSs large 
enough to account for the observed SN la rate. However, as re- 
cently pointed out, the observed number of SSSs (Di Stefano 
120 lOi) and integrated soft X-ray flux (Gilfanov & Bogdan l20T0] l 
observed from external galaxies appear to be at least one and 
more likely two orders of magnitude too low to account for the 
SN la rate. Furthermore, pre-explosion observations of the po- 
sitions of nearby (< 25 Mpc) SNe la using archival Chandra 
data have so far yielded no detections (Nielsen et al. 2012 ). This 
dearth of SSSs could very well mean that the missing SSSs are 
simply not there, and hence that the SD progenitor scenario is 
not the dominant contributors to the SN la rate. An alternative 
possibility is that the nuclear burning WDs appearing as SSS in 
SD progenitors do in fact exist and produce a significant fraction 
of the total SN la rate, but are somehow hidden from view of 
our cuiTent observational capabilities (in X-rays) during much 
of their supersoft phase. 

In the following we wish to explore the latter option. We 
consider a simple model of a massive, accreting WD with a com- 
panion star that is losing matter into the circumbinary region in 
addition to the matter it transfers to the accretor The goal has 
been to determine how much cirumbinary material is needed to 
render a nuclear burning WD in a nearby galaxy undetectable as 
a SSS for a given combination of binary parameters. 

A note on notation: in this paper we will refer to the rate of 
material that is lost into the circumbinary region simply as 'the 
mass loss rate', M. This should not be confused with the rate of 
mass that is transfeiTed to the WD accretor. Mace- In our notation, 
the total rate of mass lost from the donor is Mtot - M + M^cc 

In section|2]we describe our model, including the structure of 
the gas bubble surrounding the SSS and the contributions to the 
obscuration from neutral gas, ionized gas, and dust. In section [3] 
we present the results of our calculations, and section|4]discusses 



the observational implications of the results. Section|5]discusses 
the caveats of our model, and section|6]concludes. 

2. Model 

We consider the emergent radiation from a massive (~ IM©) 
WD accreting mass from a companion star in a close binary sys- 
tem. The donor may be a main sequence or evolved star The sys- 
tem is losing mass into the circumbinary region, and this mass 
loss has created a spherical distribution of matter (gas and pos- 
sibly dust) around the binary which may absorb and/or scatter 
the X-rays. The mechanism behind the mass loss from the donor 
into the circumbinary region is left unspecified in our study, but 
may be envisioned to be e.g. a stellar wind, wind Roche-lobe 
overflow (Mohamed & Podsiadlowski I2007I I. stellar pulsations 
of the donor, tidal interactions (e.g. Chen et al.'201 1^, a recently 
expelled common envelope, or any other process by which ma- 
terial can be deposited in the circumbinary region instead of be- 
ing accreted by the WD. It is also possible that the material in 
the circumstellar region may originate from the WD if it emits a 
wind (Nomoto et al. 1979). Whether the above-mentioned mech- 
anisms are actually capable of producing a spherical circumbi- 
nary configuration of matter is a question we do not enter into in 
this study. 

The WD is burning accreted material at its surface at the 
steady burning rate mentioned earlier. The resulting luminosity 
is that of a typical SSS, i.e. Lboi - lO"'^ erg/s, and the spectrum is 
a black body with kTu, = 50 eV. Our focus is on X-ray observa- 
tions, and as the companion is not expected to emit appreciably 
in X-rays, in observational terms our model system reduces to 
a single nuclear burning WD within a bubble of circumbinary 
material. 

2.1. The gas bubble 

We parametrize the mass loss into the circumbinary region by a 
wind velocity, m,,,, which we assume to be constant in the region 
outside of the position of the WD. We choose a value of m„, = 10 
km/s, typical for the winds of evolved intermediate mass stars 
(e.g. Panagia et al. 120061 1. 

As a first approximation, we assume the mass loss into the 
cirumbinary region to be spherically symmetric (see section |5] 
for a discussion of the caveats of this assumption). The outer 
extent of the spherical distribution of material depends on the 
wind velocity and age of the mass losing phase of the donor 
star. For the chosen wind velocity the extent of the obscuring 
gas bubble is 2.1 AUyr"' ~ lO^^* pcyr"'. 

2.2. Obscuration by neutral gas 

In general, the optical depth t along the line of sight between 
the source and the observer is the opacity k times the obscuring 
column M: 

dr - Kpdr t - kM, (1) 

where M - J pdr, and p - M /(4-Kr^Un), and M is the mass loss 
rate. 

The total neutral column along the line of sight is the sum of 
the contribution from the local gas bubble, the ISM in the host 
galaxy, the IGM between the Milky Way and the host galaxy, 
and the ISM in the Milky Way. Therefore, for a given species of 
neutral gas in our spherically symmetric model the attenuation 
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Fig. 1. Schematic drawing of the model used in this study. The 
SSS system consists of a WD accreting material from a donor 
The supersoft X-ray emission is the result of steady thermonu- 
clear burning on the surface of the WD. The WD and donor are 
surrounded by a circumbinary configuration of gas and possibly 
dust. The radiation from the nuclear burning WD may ionize a 
region around the system; in this sketch the ionized region is 
localized narrowly around the binary, but for certain configura- 
tions the ionized region may extend to or beyond the edge of the 
circumbinary gas bubble. Before reaching an observer at Earth 
photons from the SSS pass through the circumbinary material, 
the ISM in the host galaxy, the IGM, and the ISM in the Milky 
Way.|2] 



is formally given by 




Fig. 2. Opacities from K-shell ionizations for neutral gas and 
atomic metals as a function of photon energy. All values are for 
solar abundances. The dotted line gives the opacity of hydrogen 
alone, dashed is the opacity of hydrogen + helium, and solid is 
the total opacity of all chemical elements at solar abundances. 
Values are calculated from Morrison & McCammon (1983). As 
explained in section|5] the relevant photon energies for Chandra 
are from 300 eV. Furthermore, for SSS we do not expect to see 
any photons above ~ 3 keV. 



+PlSM('"host - '"1 + '"obs - ruw) 

+PIGm('"MW - '"host)jj (2) 

where is the opacity of the neutral gas, ro is the inner ra- 
dius of the neutral region of the gas species in question, ri is 
the outer radius of the spherical gas bubble, rhost is the distance 
from the source to the edge of the host galaxy, tmw is the dis- 
tance from the source to the edge of the Milky Way, robs is the 
total distance from the source to the observer, and pisM and pigm 
are the gas densities of the ISM and IGM, respectively. We put 
ro = a, where a is the separation between the binary component. 
For both the wind material, the IGM, and the ISM we assume 
solar chemical abundances. X-ray absorption happens by way 
of K-shell ionizations, and the resulting photon energy depen- 
dent cross-sections are obtained from Morisson & McCammon 
( I1983I I, as shown in figure |2] 

The formal expression in eq.© can be simplified somewhat. 
The contribution from the IGM is negligible, and we also as- 
sume that there is no significant contribution to the column from 
the host galaxy. This is routinely done for SSS studies (see e.g. 
Kahabka & van den Heuvel 2006 ), and many of the galaxies used 
in the studies of for example Gilfanov & Bogdan (2010 ) and Di 
Stefano (2010) are gas poor Therefore, the column that we con- 
sider in our model is just the sum of the column of material in 
the Milky Way and the circumbinary material. 



2.3. Obscuration by ionized gas 

Close to the source the radiation generated by nuclear burning on 
the surface of the WD will photo-ionize the hydrogen in the gas. 
Since the peak energy of a canonical SSS is low (kTht ~ 30-100 
eV) heavier elements are not likely to be appreciably ionized, 
and certainly not ionized to K-shell. 

To determine the ionization structure of hydrogen around 
the source we follow the approach first suggested by Stromgren 
(I1939I I. The Stromgren sphere is the volume around an ioniz- 
ing source in which the ionization rate equals the recombination 
rate. The rate of recombinations to atomic energy level n per unit 
volume is given by Nr^„ = nenp/3„(Te), where rig and rip are num- 
ber densities of free electrons and protons, respectively, and /?„ 
is the recombination efficiency of the n'th level, which depends 
on the electron temperature of the gas, T^. 

If we assume complete ionization within the ionized region 
then He - Up for hydrogen. Additionally, we can omit /3i from 
the expression, since every recombination directly to the ground 
level emits a photon capable of causing another ionization which 
we assume it will do immediately. The total number of recombi- 
nations per unit time is then: 

r-Ri 

r-Ri 

= 47r/?2+(r.) n](r)r^dr 
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47Tulmj 



r-Ri 

Jra 



(3) 



where j02+ ~ 2 • lO^^^r^T^^^'^cm-^/s is the total recombination rate 
of transitions above the lowest {/32+ - YjnPn - and we have 
used that 



ne(.r) = nuir) = 



M 



1 



(4) 



where nig is the mass of the hydrogen gas. 

The number of ionizing photons emitted by the SSS per unit 
time is given by: 



5, ^AnR 



Jl3.6eV 



hv 



-dv ■■ 



L 



Jl3.6eV 



hv 



'-dv 



(5) 



where h is the Planck constant, cr is the Stefan-Boltzman con- 
stant. By is the frequency dependent Planck function, and v is 
the frequency. 

Setting S , equal to Nr^ioi we get: 



(6) 



We assume that the temperature of the entire gas bubble is 
the effective temperature of the SSS. For kThb = 50eV the effec- 
tive temperature is 5.8 • lO^K. This is certainly an overestimation, 
but as discussed in section|5]the impact of this inaccuracy is neg- 
ligible. 

Clearly, eq.(l6]l only has a physically meaningful solution if 
is larger than the second term on the right-hand side. This 
gives us a constraint for the mass loss rate: 



M > 



MTe) 



= Mstr. 



(7) 



For our model system, Lboi = 10""^ erg/s and kTi,h = 50 eV, 
the number of ionizing photons is 1.46 ■ 10'*^s"', and Mstr - 

2.00-10-<'Moyr-i-(^)'^'. 

For mass loss rates larger than Mstr there is a clearly defined 
inner ionized region, outside of which only neutral matter exists. 

For mass loss rates smaller than Mstr the expression for the 
ionized sphere will not be physically meaningful. In this case, 
the assumption that all photons capable of ionizing the gas are 
absorbed and cause ionizations is inaccurate. For such mass loss 
rates the gas bubble is not dense enough to absorb all photons 
capable of ionizing the gas, and there is no longer a clearly de- 
fined ionized region. Therefore, all hydrogen and helium in the 
gas bubble is fully ionized, and excess photons seep out into the 
interstellar medium, possibly causing further ionizations there. 

In the ionized region, the only contribution to the obscuration 
from hydrogen and helium will be through Thomson scattering. 
The cross-section of Thomson scattering is largely indepent of 
photon energy, and therefore the free electrons produced in the 
ionizations will affect the absorption/scattering at all photon en- 
ergies at approximately the same amount. 

For X-ray binaries involving neutron stars or black holes 
the energies and densities involved may sometimes lead to 



Comptonization of the plasma, which has an impact on the ob- 
scuration caused by the Comptonized material. However, for the 
lower energies and densities involved in our model Compton 
scattering plays no role at all. 

2.4. Obscuration by dust 

Depending on the properties and the mass loss mechanism of the 
donor star, a fraction of the metals in the gas may be condensed 
into dust grains, and we need to consider if the presence of dust 
changes the total X-ray absorption of the circumbinary material. 

As mentioned earlier. X-ray absorption happens by interac- 
tion of photons with K-shell electrons, and in this way the ab- 
sorption cross section of an individual atom is largely indepen- 
dent of the location where the atom is found. However, putting 
atoms into grains represents a form of clumping, which can de- 
crease absorption if individual grains are akeady optically thick 
to X-rays of the considered energy. In this case, a part of the 
grain does not contribute to the X-ray opacity, because any X-ray 
photon will already be absorbed in the source-facing side of the 
grain. This can reduce the opacity by a so-called self-blanketing 
factor (Fireman 1 19741 1 of 



/,=(l-e-<V>)/T 



(8) 



Here (Tg, ) is the average optical depth of individual grains. 
However, the effectiveness of this clumping of X-ray opacity 
is very limited. Even in cold diffuse clouds in the interstellar 
medium, important elements are hardly depleted into dust grains. 
For example, nitrogen and neon are abundant elements that re- 
main entirely in the gas phase. In the mass loss flow of the donor 
star, this condensation will be even less complete. If the donor 
star has an oxygen-rich chemistry, the entire carbon content of 
the wind and an equal amount of oxygen will be trapped in the 
very stable CO molecule and not participate in the condensa- 
tion process (Gail and Sedlmayr 1986 ). For solar-system-like el- 
ement abundances (Anders and Ebihara [T982]l, this means that 
all of the carbon, all of the Nitrogen and Neon, and 40% of the 
oxygen remain in the gas phase. As these are the most important 
absorbers for soft X-rays, it is clear that the soft X-ray absorp- 
tion cross section will be reduced by a factor much less than 2. 
Wilms et al. (2000) conclude that even in the ISM where con- 
densation is rather complete, for grain sizes smaller than 0.3;um, 
the resulting effect will not exceed 10%. In the donor mass flow, 
the effect will be even smaller. Therefore, we can safely ignore 
this complication and assume that all metals are in the gas phase, 
fully contributing to the X-ray absorption. 



3. Results 

We first present the generic results of our model, i.e. results that 
do not depend on the instrument being used. Then we apply our 
results to a specific case, that of Chandra's ACIS-S detector. 

3.1. Generic results: Ionization structure and obscuration 

Figure |3] shows how the obscuration depends on orbital sepa- 
ration and mass loss rate for 350 eV photons emitted from a 
L = lO^^erg/s, kThh = 50eV SSS. Photons at 350 eV are safely 
above the photon energies at which typical X-ray observations 
are unreliable, but close enough to the peak of the black body 
curve that an unobscured SSS produces appreciable amounts of 
photons. 
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M [Mo/yrl 



Fig. 3. Contour lines showing the dependence of obscuration on 
mass loss rate and binary separation, for photons emitted at 350 
eV. The left-most solid line corresponds to an attenuation factor 
of 5, the dotted line to an attenuation of 10, dashed to an at- 
tenuation of 100, and dot-dashed to an attenuation of 1000. The 
right-most solid line is the critical mass loss rate Mstr; to the left 
of this line the wind is not dense enough to sustain a clearly de- 
fined ionized region around the source. To the right of it a clearly 
defined ionized region exists around the source, while the mate- 
rial outside of this region is neutral. 



Our study shows that for binary separations around 1 AU a 
spherically symmetric mass loss rate of ~ lO^^Moyr"' is suf- 
ficient to fully obscure the supersoft X-ray emission from our 
model system. 

It is evident from figure |3] for orbital separations around 1 
AU a clearly defined ionized region will form at considerably 
larger mass loss rates (~ lO^^Moyr"') than that needed for full 
obscuration. Hence, at this photon energy, in between those mass 
loss rates the SSS will be fully obscured, but the binary will be 
surrounded by an extended ionized region that may be detectable 
at other wavelengths than X-rays (see section|4]l. 

We note that the obscuration curves on fig|3]are for a pho- 
ton energy of 350 eV. As the photon energy rises the curves 
move further to the right in the plot, while the curve for the crit- 
ical mass loss rate remains in place (since the critical mass loss 
rate does not depend on the photon energy). This means that for 
larger photon energies one can imagine configurations for which 
these curves overlap, i.e. there is a clearly defined ionized region 
around the binary, while the X-ray emission from the binary is 
not fully obscured. However, this does not happen for our SSS 
at photon energies below ~ 5 keV, and the number of photons 
emitted by the system at these energies is miniscule. 

3.2. Specific results: Cliandra's ACIS-S detector 

As an application, figure|4]illustrates what our model system will 
look like if it resides in MlOl and is observed with the ACIS- 
S detector on the Chandra X-ray Observatory. The distance to 
MlOl is taken to be 6.4 Mpc (Shappee & Stanek 2011). The 
contribution to the column from the Milky Way is obtained from 
Dickey & Lockman 1 1990a : for MlOl we find a column of 1.15 ■ 



10^" Nn/cm^. We chose Chandra's ACIS detector as an example, 
since this is the instrument used by most groups (e.g. Voss & 
Nelema ns l2008l Roelofs et al. 2 008 Ne lemans et al. l2008l Di 
StefanolMl Gilfanov & Bogdan[2ni0l Nielsen et al.|20T2ll. 




V [Hz] 



Fig. 4. Unabsorbed and absorbed black body curves of a SSS 
in MlOl. The solid curve corresponds to a black body nuclear 
burning WD with L = 10^^ erg/s, kTi,h = 50 eV and no obscur- 
ing circumstellar matter (only galactic Nh), folded with the ef- 
fective area function of Chandra's ACIS-S detector. The dotted 
curves are the absorbed black body curves of the same source for 
seven logarithmically equidistant values of M from lO^'^Moyr"' 
to 10"^ Moyr"' (see table [TJ, with the mass loss rate increas- 
ing from left to right. The orbital separation and inner radius of 
the gas bubble is 1 .5 AU. The numbers on each curve give the 
number of photons expected to be received in ACIS-S for an in- 
tegration time of 40 ks. For comparison with figure|3]the vertical 
dashed line is at a photon energy of 350 eV. 

The photon energy dependent eff'ective area function of 
ACIS-S can be found on the Chandra homepag43. In order to 
find the number of photons detected we fold the calculated flux 
of the source with the effective area before integrating over all 
photon energies. 

Since the launch of Chandra in 1999 the sensitivity of the 
onboard detectors have degraded somewhat. We adopt the effec- 
tive area function for the earliest possible Chandra observations 
(cycle 3). If a source can be obscured sufliciently to be unob- 
servable to the ACIS detectors when these were new and at their 
most sensitive, then such a source would also be unobservable 
to the older, less sensitive ACIS detectors. We note that the de- 
tectors on Chandra are not sensitive to photons below roughly 
100 eV. Additionally, Chandra detections at photon energies be- 
tween 100 and 300 eV are known to be unreliable, and analyses 
should therefore filter out photons at energies below this thresh- 
old. This means that observations of SSS spectra are in fact only 
observations of the high energy tails of their spectra, since the 



' http://cxc.cfa.harvard.edu/cgi-bin/build_viewer.cgi7ea 
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Table 1. Attenuation (A^y,obs/M)) of the integrated number of 
photons from our model system located in MlOl, as expected 
for different mass loss rate and observed with Chandra's ACIS- 
S detector. The mass loss rates listed correspond to the dotted 
lines in figE] The integration has been done for photonenergies 
in the 0.3 to 3 keV range. For comparison, the number of pho- 
tons expected in ACIS-S from the same source integrated over 
the same energy range in a completely empty (i.e. in the absence 
of galactic or circumbinary material) is A^o = 4.4 • lO^^s"'. Note 
that the attenuation is for the integrated number of photons, and 
hence is not immediately comparable to fig|3] as the values on 
that figure are for a single photon energy (350 eV). 
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spectral peaks are located far below the detection threshold of 
Chandra. 

The number of photons we expect to receive in the relevant 
energy range (300 eV - 3 keV) for observations with ACIS-S is 
given in table[T]for each of the curves in figure|4] 

4. Observational implications 

In the preceding sections we have laid out the setup of the model 
and our results. To put our results into an observational context 
we will now review the current evidence for the existence of 
circumstellar matter around SNe la, which would be a signature 
of the SD progenitor scenario. We will also discuss the possible 
observational implications of our results for SN la progenitors. 

4.1. Evidence for CSM around SN la progenitors 

Several studies have found evidence for circumstellar matter in 
SNe la explosions: 

High-velocity features in early-time optical observations of 
SN2003du were interpreted by Gerardy et al. (2004) as evidence 
for the interaction of the outer-most layers of the SN ejecta with 
a dense circumstellar shell of solar metallicity material created 
by mass loss from the progenitor system prior to the explosion. 

In a study of SN la remnants DEM L238 and DEM L249 
Borkowski et al. (2006) found bright central X-ray emissions 
surrounded by fainter shells and interpreted this as remnants 
of circumstellar media around the progenitors that had been 
shocked to emission by the SN ejecta. 

In a more recent study, Chiotellis et al. (1201 11 1 compared 
2D model simulations with observations of the historical SN la 
SN1604, also known as Kepler's SN. The SN remnant shows 
a peculiar nitrogen-rich shell-like structure in optical images. 
Simulations by Chiotellis and collaborators assumed a mass loss 
of 10"^ - lO-''M0yr-' and wind speeds of 5-20 km/s, typical of 
thermally pulsating asymptotic giant stars. The observed shell- 
like features are reproduced by their simulations, and interpreted 
as a shocked interaction layer between the progenitor's wind- 
blown circumstellar bubble and the SN ejecta. If correct, their re- 



sults are consistent with a SD progenitor emitting a stellar wind 
prior to explosion. 

Another recent study by Patat et al. (120071 ) found blue-shifted 
absorption features of the Na I doublet (5896 A and 5890 A) in 
optical spectra of SN2006X, which was interpreted as evidence 
for gas outflows from the progenitor system. Their results were 
generalized by Sternberg et al. (1201 It who found similar features 
in an unbiased sample of 35 SNe la, indicating that such features, 
while not demonstrably present for other types of SNe, may be 
characteristic of SNe la. 

4.2. Upper iimits 

Conversely, a number of studies have looked for, and failed 
to find, evidence for circumstellar matter in SN la explosions. 
These non-detections have led to upper limits being placed on 
the possible mass loss from SN la progenitor systems in several 
wave-bands. 

Using optical spectra of SN2001el, Mattila et al. ( 120041 ). 
found upper limits of M < 9 ■ lO"*^ and 5 ■ lO-^Moyr-' for the 
progenitor system of the SN, for wind speeds of 10 and 50 km/s, 
respectively. 

Studies aimed at finding radio emission from the interac- 
tion between SN ejecta and circumstellar matter have been un- 
dertaken by several groups. No direct detection has been made 
at this point, but upper limits are reaching interesting values. 
Using the VLA, Panagia et al. 2006 found upper limits of ~ 
lO-''M0yr-' based on observations of 27 SNe la. Their analysis 
extrapolates from the assumption that the process behind radio 
emissions from SNe la are similar to those of SNe Ib/c. 

Using more recent radio observations, Chomiuk et al. ( 1201 II) 
analyzed EVLA observations of early SN la spectra. From 
their non-detections they found typical upper limits M/u^. < 
lO-'M0yr-V(lOOkm/s) for most sources (private communica- 
tion). In another recent article, Chomiuk et al. (2012) reports up- 
per limits of M/u^ = (6 • 10-'" - 3 ■ lO-'')M0/yr/(lOOkm/s) for 
non-detections of radio emission from SN201 Ife, the closest SN 
la in 25 years. For the lower wind speed used in our model the 
upper limits are correspondingly smaller; hence, the upper limit 
on the wind mass loss rate of the donor becomes lO-'^M0yr-' for 
typical SD SN la progenitors, and (6 ■ 10-" - 3 ■ lO-'")M0/yr 
forSN2011fe. 

4.3. Observationai predictions 

We note that except for the case of SN2011fe the upper limits 
found in the studies mentioned above are all larger than what we 
require for full obscuration of systems with binary separations 
of < lOAU. The upper limits found by Chomiuk and collabora- 
tors come closest to constraining our results, and if their general 
limits are correct then our model cannot explain obscuration of 
systems with binary separations larger than 1 AU. If their limits 
for SN2011fe are correct then the non-detection of X-ray emis- 
sions from that particular SN cannot be explained by obscuration 
from circumbinary material in our model for a giant donor, since 
that would require binary separations of < 10-^ AU, effectively 
placing the WD within the envelope of the giant. A MS donor 
cannot be ruled out by these radio upper limits, neither can a 
WD wind. Hopefully, future observations will provide either a 
detection of the shocked region or stronger general constraints 
with which to compare our model. 

Another important point is that the density (and thereby mass 
loss rate) is not the only important parameter when determin- 
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ing upper limits. Our analysis shows that even relatively small 
circumbinary gas bubbles are able to fully obscure the system: 
As can be seen from the ro-dependence in eq.Q, the obscura- 
tion is caused mainly by material in a very narrow region around 
the SSS. The outer extent of the circumstellar gas bubble is es- 
sentially irrelevant, except for very small bubbles where ro is 
comparable in size to ri. This means that even quite compact 
systems, say up to a gas bubble radius of ~ 10 AU, would be 
capable of obscuring a SSS, provided the mass loss rate is large 
enough. It follows that the mass loss does not have to have been 
'on' for very long to obscure the system in X-rays. We empha- 
size that this does not solve the problems raised by Di Stefano 
(120 lot and Gilfanov & Bogdan (2 010) . as the source still need to 
be 'on' for a significant period (~ 10^ years) to be able to grow 
significantly in mass. But it may explain why many SSS appear 
to be highly variable. 

In general, observations to detect or constrain the mass loss 
rate of the progenitor systems of SNe la need to be performed 
a very short time after the SN explosion. The bulk of the SN 
ejecta moves at ~ 10, 000 km/s, corresponding to roughly 6 AU 
per day. Therefore, the interaction shock from a nuclear burning 
WD shrouded in an fully obscuring gas bubble with a radius in 
the range of a couple of tens of AU is unlikely to be detected by 
anything but the very earliest (1-2 days after explosion) observa- 
tions. Since radio appears able to supply the best upper limits the 
ideal observing scheme would be to obtain EVLA observations 
of SNe la within a day or less of the explosion. 

If the companion in our model system is evolved we can dis- 
regard orbital separations less than ~ 0.5 - lAU, since that will 
be inside the outer layers of the companion. At such short sep- 
arations the WD is more likely to spiral into the companion or 
cause the envelope to be expelled, rather than go through a stable 
period as a SSS. However, several studies argue against the pos- 
sibility of giant companions. For example, using pre-explosion 
archival Hubble Space Telescope (HST) images Li et al. ( 1201 11 1 
ruled out a luminous giant or supergiant as the companion to 
SN2011fe, although a sub-giant companion could not be ex- 
cluded by the data. 

If the companion is a main sequence (MS) star the system 
can exist stably at much smaller orbital separations, and the mass 
loss rates required for obscuration are correspondingly lower, as 
evident from figure [3] But the mass loss rate from such systems 
may also be much lower than for systems containing evolved 
stars. 

There is also the possibility that the mass loss from the sys- 
tem is caused by a wind from the WD itself. This could happen 
if the mass loss from the donor to the accretor is slightly larger 
than the maximal steady burning accretion rate. The accreted 
mass could be supplied by either a MS or evolved donor If the 
accretion rate is larger than the steady burning rate, the WD will 
'puff up' from the accretion and emit a wind of its own, while 
possibly still burning some of the material on its surface, see 
Nomoto et al. ( Il979l l. The 'orbital separation' in figure [3] should 
then instead be perceived as the diff'erence in radius between the 
nuclear burning layer and the wind emitting layer of the WD. 
For such small separations full obscuration can be achieved even 
with fairly small (10"'' - lO^'^Moyr"') mass loss rates. If we 
envision a WD accreting at slightly above the maximal steady 
burning rate and emitting a weak, spherical wind of this magni- 
tude such a source would be completely obscured in our model. 
Potentially of more general interest, if it could be shown that the 
accretion process feeding a steady burning WD is not 100% ef- 
ficient, or if the X-ray emitting surface of the WD loses a very 
small fraction of the accreted material to the circumstellar region 



while it burns the result could be significant or complete obscu- 
ration. Therefore, these WD wind scenarios could potentially ex- 
plain the absence of a large fraction of the SSS that become type 
la SNe in the SD scenario. 

A possible observational characteristic of our model could 
be H-ff emission caused by recombination in the ionized gas 
bubble. This emission may be visible in archival optical images 
of nearby SNe la, and it would be logical to suggest a system- 
atic archival search for this kind of emission in pre-explosion 
images at the position of nearby SN la progenitors. However, 
such a search of the HST archive was performed by Voss et al. 
(2012, in prep.), who found no evidence for optical counterparts 
for nearby SNe la progenitors, so it is unclear whether such a 
search is actually feasible. Another option is to search for H-a 
emission in regions where such emissions are not to be expected, 
i.e. outside of young, star-forming regions. It might be possible 
to make an analysis similar to the SSS 'counting' done by Di 
Stefano (2010), but in H-a instead of X-rays. We note, however, 
that even though dust is unimportant for the obscuration of X- 
rays from our model system it could possibly obscure the H-a 
emissions, since dust is a more efficient absorber at optical wave- 
lengths. Depending on the tenuousness of the outer parts of the 
gas bubble there may also be forbidden lines, but this depends 
heavily on the density of the gas. 

4.4. Symbiotics 

Our results may also explain why few classic symbiotic binary 
systems are not observed as SSS. Despite the fact that the WDs 
in symbiotic systems are expected to be massive and accret- 
ing mass at rates comparable to the steady burning region only 
three symbiotic SSS are currently known (SMC3, Lin 358 and 
AG Dra). Somewhat analogous to the evidence for circumbi- 
nary matter around SN la progenitors observations of outbursts 
from recurrent novae in symbiotic systems also indicate the pres- 
ence of significant amounts of circumbinary material, provid- 
ing absorbing columns large enough to fully obscure the sys- 
tems in quiescence (e.g. Shore et al. 1996). The physical param- 
eters of typical symbiotic systems are Lboi - 10^ - 10*Lq = 
4 ■ 10^^ - 4 ■ lO^*" erg/s and a = 2 - 5 AU, while the red giant wind 
emitted by the donor in such systems is roughly 10"^Mo/yr, see 
e.g. Mikolajewska (I2012I I. If we assume that the spectra of this 
type of source is a black body comparable to a canonical SSS 
then Figure [3] shows that even if the luminoisities of these sys- 
tems were as large as the ones expected for the nuclear burning 
SD SN la progenitors they would be completely obscured. The 
fact that they are observed to be one to two orders of magnitude 
less luminous only serves to make them even easier to obscure. 

5. Discussion 

As mentioned earlier, we use a number of simplifying assump- 
tions in our calculations. Here we discuss the caveats introduced 
by these assumptions. 

5.1. Density profile 

The assumption of constant wind speed from the surface of the 
companion star is probably not correct. In reality, the wind is 
accelerated by a variety of processes until it reaches its terminal 
velocity, and this is not expected to happen until well beyond the 
orbit of the binary. Therefore, our wind speed is probably too 
large. If the wind is accelerated through the system and doesn't 
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reach the constant value used in our simulations until some time 
later, i.e. further away from the source, the density close to the 
source will be larger. So, this assumption also underestimates the 
amount of obscuration. 

Also, we have assumed a spherically symmetric distribu- 
tion of matter around the binary. This is in general not what 
is observed in symbiotic systems, where disk-like structures in 
the orbital plane are expected, possibly accompanied by bipo- 
lar outflows (e. g. Solf &Ulrich [T985t Corradi & Schwarz [T993l 
Munari & Patat'1993i Overall, it is difficult to say whether the 
assumption of sphericity is likely to over- or underestimate the 
amount of obscuration. In case of a non-spherical structure, the 
obscuration in specific cases depends sensitively on the incli- 
nation to the sight-line of the observer. This uncertainty could 
be dealt with if we had a firm understanding of the symmetries 
of the matter in the relevant binary systems. Given the absence 
of that understanding, combined with the earlier mentioned fact 
that only the material very close to the system has a significant 
effect on the obscuration, we believe that our spherically sym- 
metric model is a reasonable first approximation. Hopefully, fur- 
ther studies will provide a better understanding of the density 
structures of the systems in question. 

5.2. Temperature 

In our calculations we have a assumed a constant temperature 
throughout the entire gas bubble, i.e. the temperature of the sur- 
face of the WD. In reality, the temperature of the gas bubble will 
fall off with the distance from the emitting SSS. However, since 
even at a temperature corresponding to the Thi, of the SSS the 
elements of importance to absorption and scattering in the ob- 
servationally relevant interval, i.e. elements heavier than helium 
but lighter than iron, will not be fully ionized. Our simplified 
temperature assumption only plays a role for hydrogen and he- 
Uum, which are unimportant absorbers at the photon energies 
where Chandra is sensitive. We therefore estimate this effect to 
be negUgible for our purposes. 

5.3. Dust and stellar winds 

As explained in section l24l dust appears to be fairly unimportant 
in connection with the possible obscuration of a SSS. 

However, if the X-ray source manages to ionize a large re- 
gion in the circumbinary gas, in this region dust formation will 
not be possible. If radiation pressure on dust is an important fac- 
tor for driving the wind in the first place, introducing the X-ray 
source may have significant effects on the mass loss rate from 
the companion. If dust has formed already, it might be destroyed 
again, for example by thermal evaporation due to X-ray heat- 
ing of the grains (e.g. Fruchter et al. 1200 1 1) , by charge explo- 
sions due to massive photo-electric ionizations (ibid), or by ther- 
mal sputtering in hot gas (e.g. Tielens et al. l200Il l. In this study 
we have not considered such effects in detail. Instead, we have 
used a fixed mass loss rate as a parameter of the model and dis- 
regarded the consequences for X-ray absorption if part of the 
heavier atoms are present in the form of dust grains. 

5.4. Metallicity 

By using the model of Morrisson & McCammon (1983) we as- 
sumed solar metallicity of the obscuring material (both the wind 
material and ISM). As mentioned, the most important obscur- 
ing elements for the photon energies accesible with Chandra are 



elements heavier than helium but lighter than iron. For metallici- 
ties different from the one used in this study the obscuration will 
scale accordingly. 

5.5. Spectrum 

Following the original work of van den Heuvel et al. 119921 we 
have assumed the SSS to be a simple black body. This is prob- 
ably not an entirely accurate description of a nuclear burning 
WD (e.g. Ness et al. 2003 Rauch 2003 Rauch & Werner|20T0ll, 
and by assuming a black body spectrum we may well overesti- 
mate the temperature of the actual surface of the WD. However, 
deviations of the actual spectra from that of a black body is neg- 
ligible in this context, since we are analyzing obscurations of 
several orders of magnitude. In addition to this, the temperature 
dependence in our calculations is quite weak (e.g. the number 
of ionizing photons in the calculation of the extent of the ion- 
ized region depends on temperature like ~ T"^^"*). Observations 
of SSS have often assumed a black body spectrum (e.g. Greiner 
2000), so comparisons of such observations with our model will 
in a sense be consistent. We note that the question of WD at- 
mospheres is not particularly well understood at this point, so 
a very detailed analysis with additional assumptions would not 
necessarily improve the applicability of our results. 



6. Conclusions 

To date, it has mostly been assumed that nuclear burning WDs 
in SD progenitor systems would be more or less 'naked'. 
Consequently, the absence of a large enough number of these 
sources have been seen as a problem for the SD model, as there 
seem to be too few of these sources to account for the observed 
SN la rate. 

We have examined a model system of a canonical SSS em- 
bedded in a spherical circumbinary gas bubble. The mechanism 
behind the formation of the gas bubble has been left unspecified, 
but could be the result of e.g. a stellar wind from an evolved 
companion, wind-RLOF, pulsations of the donor, or tidal effects 
between the binary components. 

We have shown that for a certain critical mass loss rate (e.g. 
M ~ lO^^Moyr"' for a ~ 1 AU) a clearly defined, narrowly situ- 
ated, ionized region will form around the SSS. For systems with 
mass loss rates below this critical value the SSS will be sur- 
rounded by extended ionized regions that may extend into the 
ISM. 

Our results suggest that for systems with a ~ lAU quite 
modest circumbinary mass loss rates (~ 10"'' - lO^^Moyr"') are 
sufficient to significantly obscure the nuclear burning WD SSS. 
This mass loss is in addition to the mass the donor loses to the ac- 
creting WD. For wider systems larger mass loss rates are needed. 
Even at orbital separations on the order of 100 AU the mass loss 
rates required for significant obscuration (~ lO^^Moyr"') are not 
unrealistic for some late red giants or asymptotic giant branch 
stars. However, if they exist such wide systems are unlikely to 
be able to supply the mass loss rate required for steady burning. 

The mass loss rate required for the ionized region to be- 
come clearly defined is several orders of magnitude larger than 
the mass loss rates needed for total obscuration at the relevant 
photon energies (< IkeV). This means that SSS systems with 
sufficient mass loss rates to fully obscure the X-ray emission 
will have an extended region of ionized material surrounding it. 
According to our model, for binary separations ~ lAU mass 
loss rates between 10"^ and lO^^Moyr"' produce such systems. 
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While not observable as SSS, these systems may instead be ob- 
servable in IR, radio or H-a as recombination nebulae. To the ex- 
tent that they are available multi- wavelength archival searches of 
pre-explosion images at the positions of nearby SNe la may find 
such emissions, even if no sources have been found in archival 
Chandra images. We also propose that H-a emissions outside of 
young regions may be evidence of obscured SSS. 

If a steady burning WD emits a small amount of the accreted 
material, on the order of ~ 10"' ' - lO^'^'Moyr"', our calcula- 
tions show that they may be completely undetectable in X-rays. 
We have no way to determine whether steady burning WDs ac- 
crete at 100% efficiency, but if it does not then even such a small 
amount of material will have important consequences for the X- 
ray signature of such objects. 

The full obscuration constraints for binary separation and 
mass loss rate presented above are probably too strict. Our model 
examined an observational best-case scenario, and when simpli- 
fying assumptions have been made they have consistently been 
made to favor a minimal amount of obscuration. For these rea- 
sons, even lower mass loss rates may be sufficient to fully ob- 
scure more realistic systems. 

Our results may have implications for the SD scenario for 
type la SNe. The fact that it is comparatively easy to hide the 
X-ray emissions from nuclear burning WDs may help to explain 
some of the 'missing' systems mentioned in ffie Introduction. 
If a significant fraction of ffie progenitor systems can be shown 
to be severely or completely obscured by circumbinary mate- 
rial originating in the binaries themselves ffien it would explain 
the discrepancy between the SN la rate and the low number of 
observed SSS systems and integrated X-ray luminosity of ellip- 
ticals. 

Our study may also explain why so few symbiotic systems 
are visible as SSS, since typical symbiotic systems are embedded 
in a dense wind, and are less luminous ffian ffie expected SD SSS 
systems. 

In future work we plan to include our model in population 
synthesis simulations to determine if systems with the parame- 
ters required for obscuration are produced in large enough num- 
bers to account for a significant fraction of the 'missing' SSSs. 
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